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INTRODUCTION 


I i 

Galactic nebulae which show emission-line spectrum are of two 
types, that is, diffuse and planetary nebulae. The diffuse nebulae are 
normally irregular, often of low density and surface brightness, and 
sometimes extensive; whereas, planetary nebulae are generally symme- 
trical and compact and usually have a higher surface brightness and 
density than the diffuse bright nebulae. The spectra of these gaseous 
nebulae as well as those of certain stars, including the Sun, show 
strong emission lines owing to allowed and forbidden transitions of a 
number of "heavy’ 1 elements (atoms and ions other than H and He) , such 
as N, 0, Ne» S, Cl, Ar, P, Fe, Ca, Mn, Cr, V, Ca and Ni in various 
stages of ionization. In addition Mg I lines have been observed and Na 
and K are known to exist in planetary nebulae although not in great 
abundance . ^ 

Interpretation of gaseous nebulae requires proper evaluation 
of their geometrical structure , physical state, and chemical composition. 
One must be able to interpret the spectra quantitatively in terms of 
physical parameters and chemical composition of the emitting medium. 
Consequently, collisional cross-sections for allowed and forbidden tran- 
sitions and line intensities are necessary for the proper interpretation 

of the properties of gaseous nebulae. One type of excitation mechanism 

* 

which is possible ise + A-+e + A. 

This paper considers excitations of neutral atoms by inelastic 
scattering of incident electrons. Simple antisymmetrlzed product Slater 



Wave functions describing the initial and final states of the atom are 
used. Total cross-sections using the Born Approximation are calculated 
for the following excitations: Li(2s 2p) , Li(2s 3p), Na(3s 4p) , 

Mg (3p 4s), Ca(4s -*• 4p) and K (4s + 4p) . Since intensities of the 

. emitted radiation are also of importance in interpreting gaseous nebu- 
lae, the line strength of the Na(3s 3p) transition is calculated using 
Slater Wave functions and compared to an experimental result. 

4 

Typical electron temperatures in gaseous nebulae are 10 °K. 

If a Maxwell distribution is assumed the average kinetic energy at this 
temperature is about .8 ev. This is below the threshold energy of all 
of the transitions calculated here. There are of course electrons with 
energies higher than this average and these appear in the tail of the 
velocity distribution. In particular for a Maxwell distribution at 
10^ °K, .9% of the total number of electrons have energies higher than 
5 ev. This is above the threshold of all the transitions considered 
here, the largest threshold being 3.75 ev. 

Due to the relatively low abundance in planetary nebulae of 

I 

the elements for which the total cross sections are calculated here, the 
particular cross sections given may not be of great interest to those 
studying these nebulae. However, by using either of the programs given 
in the Appendix it is relatively simple to obtain cross sections for 
any desired s -*• p transition. All one needs are the threshold energies 
and the configuration of the atom to allow the determination of the 
effective atomic numbers. These effective atomic numbers . are calculated 


J 





using Slater’s rules which are given in the next section. Any deter- 
mination of the effective atomic numbers using these rules takes only 
a few minutes. 

The reason for choosing Li, Na and K was so that a compari- 
son could be made with another readily available calculation, i.e. 

2 

(Vainshtein, Opykhtin, Presnyakov). Ca and Hg were chosen because 
they were originally thought to be of importance in the emission spec- 
tra of comets. 

The method is not expected to yield reliable results near 

l ' 

threshold but should give answers within a factor of 2-7 for incident 
electron energies about 3-4 times threshold. 
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WAVE FUNCTIONS 

Two types of wave functions are used. Both have the general 

form ip . = R „ (r)Y„ (ft) where R n (r) is a Slater radial wave function 

r n£m nZ ' £m nil 

and Y^(ft) is a spherical harmonic. In one type which will be referred 
to as Slater Wave Functions I (SWFI), the ^.-dependence of the R^(r) T s 
was deleted. This is equivalent to setting l- 0 in the radial func- 
tion for all n. The general form of the SWFI is given by the following: 





i 


To make the radial integrations simpler the radial wave func* 
tions may be rewritten in the following form, 


R (r) = N l F , r n_j exp(- “ ) 

n n V , n , n+l-i ^ n+lH a 

3=1 J J o 


( 1 ) 


where the F .'s can be determined from 
n>j 

n+l-j 

F = N T F A 

n,j "j L ml n,n+l-X n+l-K,j 


( 2 ) 


and we choose F = 1, A = 0. (2 1 ) 

n,n n,n 

The Z ’s are effective atomic numbers which are determined by 
n 

4 

the following Slater rules. 

"For determining » Z-s , the electrons are divided into the 
following groups, each having a different shielding constant: Is; 2s, p; 

3s, p; .3d; 4s, p; 4d; etc. (The s and p are grouped together but the d 
and f are separated.) The shells are considered to be arranged from 
the inside out in the order named. . The shielding constant s for any 
group is found as follows, 

1. Nothing from any shell outside the one considered. 

2. An amount 0.35 from each other electron in the group consi 
dered (except in the Is group, where 0.30 is used instead). 

3. If the shell considered is an s,p shell an amount 0.85 
for each electron with total quantum number less by one, and an amount 
1.00 from each electron still further in; but if the shell is a d or f , 
an amount 1.00 from every electron inside it." 



The second type of wave functions which are used allow for an 


^-dependence in the radial function. This is accomplished by writing 
different radial functions with different ^-values for the same n. These 
radial wave functions are referred to as Slater Wave Function IX (SWF II) 
and are given by the following: 

R 10 = N 10 exp (_Z l r/a o ) 

R 20 " N 20 l- r exp ^ -Z 2 r/,2a o ) + A 201 R 10 ^ 

R 21 * N 21 & exp <- Z 2 r/2a o^ . 

. n-1 

R n* = N n£ b n_1 exp (-Z n r/na Q ) + I A , R ", ^ ] 

n =£+l 


The A n ^ n ! r s are the amount of the R n , *s with n 1 < n - 1 which must be 

added to R to make the R 's orthonormal. The N ‘s and the A . ’s 

n£ n£ n£n 

are again determined from the orthonormality conditions 


<R 6 , F f > = 6 , 

n£* n l nn’ 


n+n’+l 


A ^n' =_N -'° &n+n'): 


n £ 


Z , Z 
n , n 
~ + ~ 


n'-l 

“ £ A »tK.Vj” >n 


2n+l 


and 


N 


n£ 


[ (2n) .* 


2Z 


n 


n-1 0 u 

i ‘LjiC 

j=l J 


n 


Here again the radial functions can be put into a form which make the 
radial integrations simpler. Thus. „ 
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BORN APPROXIMATION 


The Schrodinger equation for the atom plus incident electron 


is 


2 2 2 Z 2 Z 2 

|_ v + ^_ l v.2 + E + Z^_ + J 

2m r 2m , - i T r 4 , r 
i=l J=1 j 


- 1 


l r i - r 1 


2 

-I 


3 =1 r - r ; 


¥ (r,r ) - 0 
6 


(4) 


where is the total energy, 

A 

r refers to the position of the incident electron, 

■ y 

r, refers to the position of the incident atomic electrons, 

1 

r = (r, . . . r 0 >. 
e 1 Z 

-y -y 

Now expand ^Cr,^) in terms of a complete set of orthonormal 
atomic wave functions ip^(r^) satisfying the equation 


*2 Z 

2m l 


*1 


2 

Ze Z 


1=0 


j =1 r j 


~hi 


= + E 

i-> -> n 

r. - r. 


fci < r J 

n e 


= 0 (5) 


where E are the associated eigenenergies. Now, set 


1'(?.? e ) = S t r t iS) <t A (r e ) (6) 

where S denotes a summation over the discrete states and. an integration 
over the continuum states. 

* 

■ 

Substitute equation (6) into equation (4), multiply by 'iff (r ) 
and integrate with respect to r^ to obtain 
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^ v 2 + | if v . 2 + B T+ ^ + ! 


Z 2 Z 2 

f Ze w y e 

4 - r. v, . i*> ■* 

3 M*3 l r i " r j 


V 2 

" I ®_ 

j=l i ** 

J r - r 


S/^(?) 'l J £ ^ e ) x t n ^e ) d *e = ° (7) 


Now use equation ( 5 ) to obtain 


%- [ V 2 + K 2; ] F (?) = 

2ra L * r n J n 


■f< 


(r ) V (r,r ) ¥(r,r ) dr (8) 
e € e e 


wnere 07 9 

lr ,* + , Ze 2 , r • e O) 

V(r »r_) = - — + l — __ 

e r j=i 

is the interaction energy between the incident electron and an atom com- 
posed of z electrons and- protons, and the wave number K is given by 

K » 2 = If ’ «I - V 

•'ll 

*>■ j « « 

Now suppose the incident electron has wave vector and impxnges 

“h 

on an atom in state with energy E. . ‘ Then the asymptotic behavior of F f (r) 

X ^ 


for larger r has the form 


F f (r’) * e 


iK i r< 


iK r r 1 
f 


S fl + r’ 


f f (e.O 


where f f (0,<fO is the scattering amplitude corresponding to the excitation 
of the fth state of the atom and 9 ,<J> are the polar angles of r referred 
to the direction of incidence as polar axis. 

Consider the .equation for the Green’s function 


£ V 2 + K G (?,?’) = 6(? - ?’) 
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then the solution, to equation (8) can be written as 


F f (r*) - e 


iK. * r * 0 

1 «-. + % 
fl * 2 


G(r,?') / f (r ) V (?,r a ) V (r,? e ) d^dr (11) 


In order that F should have the correct asymptotic form given by equa- 
tion (10) we find^ 


1 + 


iK„ r-r 


G( r » r ) “ 4 7T | 1 . 


( 12 ) 


where K is a wave vector in the direction of r. 

£ 

Now 


| r-r’ | = (r 2 + r' 2 - 2rr’ cos 6)** = r 1 [(-p) + 1 - 


2r Q “} 

r T cos 6J . 


So for large r’ 


|F^"¥'| * ~r 


(13) 


and 


K ' |? - P | MC r' [l - cos 0] = K f rV - K f -r (14) 


Then 


iK.r ' -± 

_ f -iK -r 

_ , + _ e e f 

G (r,r ) ^ ; 

47tr 1 


(15) 


Substitute equation (15) into Equation (11) to obtain 


it .P 


F f (r') = e 


rr — iK *"r 

«fi - 2^2 ^ h f *f (? e> v < 16 > 


m 


Comparing equation (16) with equation (10), we see that 


m 


f f ( 0 .o = - 2^2 


rr -iK 

e 4> r (r ) v (r,r ) ^ (r,r ) dr dr 

Y f e e e e 


JJ 


(17) 





To obtain the Born scattering amplitude make the assumption 
that the incident electron interacts only slightly with the target 
atom so that its wave function may be closely approximated by the 

i K T ■» 

plane wave e *■ which would be the correct wave function in the 
absence of all interaction. This approximation should be valid when 
the speed of the electron is great in comparison to that of the elec- 
tron in the target atom. This is essentially equivalent to the re- 
quirement that K^a>>l where a is the range of the electron-atom 
interaction. With the above . assumption substitute 

iK.r 

! ' T<?,r e ) = e 1 *.(?,) 

into equation (170 Then 

f £ (8,*> ■= - — f— f e i<S i - V fi (?)d? (18) 

2ith j 


where 




(? e ) V (?,? a ) ^(? e ) d? e 


(19) 


Now substitute equation (9) and equation (19) into equation (18) and use 
Bethe’s integral^ 


] 


I 


-+■ -*■ 
iK-r 2 

i ,**• Att 

* ~ dr 2 = ~2 
: r l" r 2 K 


iK’r. 


to obtain 



a 


-12 


f ^ (e ,4») = 


K 


* v * / iK*^ 

*£ *i (r e } V e / 


2m. 


dr 


2 2 
K -tt 


j=l 


/ *; < j , 


•*V 

iK*r , 


* -f 


) e ■ dr 


where K - and dik is the momentum change of the incident electron. 

The first term vanishes due to the orthogonality of the atomic 
wave function. 


So ( in atomic units) 


f f ( e ,$) ~~2 


K 


L I 


i K # it 

rp (r ) e ^ i|».(r ) dr (20). 


■ y m y 

Equation (20) can be simplified by writing ip. (r ) and \p (r ) in terms 

I C 1 s 

of antisymmetrized product wave functions. Thus 

+f<*.> ■ Ar-A = (z:)_!s e ij... K 


and 


♦Ac 1 = WA - <Z:r!S £ z m ...n *A> +A>-*n<V 


where e . . is the permutation symbol 

1J i • « K 


+1 even permutation 


e._, v - \ 0 any two indices the same 

* # » K ) 




-1 odd permutation 
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Consider the expression 


->■ 

* - iK ' r t * + 

* f < r e > e W dr e 


-> -+ 
iK*r 


dr z * (r r 


?) e C <C i (? 1 --.? z ) 


I e ij . . .K £ to...n ( d? r-- d? Z elK ‘ rt ^ ( ? 1 ) --^ Il ( ? z ) (21) 


In order to simplify equation (21) consider the case t = 1. 


Then equation (21) becomes 


( -+■ -+■ 

r * + lK.r- 

” — . e.. „ e. |dr_ ip. (r. ) e ip n (r,) 6 ...6 V 

Z ! ij . . . K &m...n | 1 1 % 1 jm Kn 


where $.(r.) is the final state wave function that does not appear in the 
ll 


initial state wave functions. 


-T 

r ->■ * -y r i 

j... K J dr ih .< r i> e V r i } <22) 


15 Z! e ij...K e £j 


Now use the relation e A e f „ = (P-1) I 6 f where P Is the 

a$y...7r a'3Y...7T aa' 


number of subscripts on the c „ 's« 

Ot PY • • * TT 

Then equation (22) becomes 


= 1- < z 


f * -v 

_1):6 u J dr iV < r i>. e \ (r i ) 

f . * - iK - r l , 

J dr ih <v e 'f i (r P 
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But thi's shows that all of the terms in the sum in equation (21) are 
equal. And since there are Z terms in the sum equation (21) becomes 


* ih* r 


dr ^ (r) e 


(r) 


where ^ (r) and \| ; ^( r ) are the final and initial state wave functions of 
the excited electron. 

Using Slater- Wave functions the expression for the Born scat- 
tering amplitude becomes 

1 /'• —y 

fO^) = - K dr'r 2 R (r) R (r)<£ f m f |e lK ’ r ! 

K J f f 2. 1 

where 

^ —V f* # Hj: 

<£ r m r |e lK ’ r ) = dft Y* (fi) e lK ’ r Y (fi) 

' f f ' ' x l f i»m. 

J f f 

To do the angular integral choose K to be along the ?-axis and 
use the expansion 

->*><» 


e ^' r . I i V j v (Kr)Y vo (fi) 

v=0 


where j (Kr) is a spherical Bessel function of the first kind. 


For an s -*■ p excitation the angular integral becomes 

= i /3 (Kr) 


<1 m- I e 
f 


iK* r 


oo> 



To do the radial integral for SOT I use equation (1 ) which gives the 

R 's in the terms of the F's. 
n 

Finally the scattering amplitude for an s -> p excitation 
becomes (for SOT X) 


2/3 N N 
n r n. 

f oo l 


. F 


(a jK~ 2)I 


4* 

a i 


. n^n +l~j n. ,n.+l-K a 

3“1* • «n f f f 11 p 

K*l. . .n^ jK 


sin (a jK 6 jK -TT/2) + cos [(a jK -l) 6 jK -*/2l]j (22) 


where 


a.-» n- + n. + 2 - j - K, 
3 K f x 


(23) 


2 n OT 1 -j .^n.+l-K 

b = [ + -a ] , ( 24 ) 

3 - • n.+l-K 

n f +l-j 1 


0 = tan 1 ( t ~~~ — ) and p = J b 2 + K 2 

jK jh jK 


For an s -> p excitation the angular integral remains the same when SWF IX 

are used. In this case use equation (3 ) for R A (r) and R . (r) . 

n^* - n . x. . 

ff xx 


So for an s -*■ p excitation with SOT II 
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‘ 2 /3N N 

n.i6. 

f(K) - — — T E • . , . F „ , -I ir 

4. j=l, . "fVY 1 * 51 VW 1- " 

K x ft 

K=l. . .n.-£. 

x x 


j a". ^ Sin (a iK 6 -iK‘ ,r/2) + p iK cos C (a iK~ 1)0 iK' ,r/ fl } (25) 


jK 


L 


jK jK 


jK 


jK ^ jK 


jK 


Total cross sections were also calculated using an effective 
potential of the form 


V ( r - r* ) = 


In dimensionless units 


V (|? - r'|) - - 


-C K 2 j r-r 1 


i-*- 

r-r 


I *')■ ”4“ 

1 *3 2 r-r 1 

-C(- y { Ka )^( | r ^ r 1 ) a n 
a o a o 

o o 


r-r 


V ( ? - P ) 


-(f-) K' 2 (I?-?' I)' 


*> *>• 

r-r 


Now choose 


c _ 7 2 „2 

a ~ Z E T 
o 


where Z is the atomic number of the target atom and E is numerically 
equal to the threshold energy in Rydbergs. The effective potential has 
an overall effect of replacing the Born scattering amplitude by 

f(K) = f Born (K)/CL + z 2 e 2 k 2 ) 



TOTAL CROSS SECTION 


The total cross-section is given by 


K r / . 

o T = |f ( e ,<<0 | dfl 


Since our V(r) has azimuthal symmetry 


dQ = 2ir sin 0d9 


the total cross section becomes 

2it A 

a T K 


| f (0) | 2 sin 0d0 


Now from conservation of momentum (hK. is momentum vector) we get 


K 2 « K. 2 + K 2 - 2K.K- cos 6 
x f if 


where hK is the momentum transfer and 8 is the a ng^^j between K_^ and 


■f 

Then 


or 


Now at 8 =0 


2KdK - 2K.K, sin 6d0 
i f 


* KdK 

sm Od0 = '-—r 
K « 

X f 


(K . ) 2 = (K. - K ) 2 or K . - K.-K 

min i f min x f 



(K maX ) 2 = (K. + KJ 2 or K maX = K + K . 

I X II 

From E = K 2 and E_ = E, - where E* is the energy of the scattered elec- 
f l T f 

tron, is the incident electron energy and is the threshold energy 
we get 


v max 

K 

mm 



± K, 





[i ± a - 


¥7 



And so the integral for the total cross section becomes 

K 



min 


(26) 


Slaters 7 rules are used to determine the effective atomic num- 
bers 2 's. These values are used to calculate the orthonormality constants 
n 

A 1 s and . N 7 s for the SWF I and the A 1 s and N 's for SWF II. The 
nj n n£n n& 

F 1 s and F f s are in turn calculated from these orthonormality constants, 
nil n£j , 

These results are then substituted into equation (22) for SWF I and into 
equation (25) for SWF 11. The scattering amplitude thus obtained must 
be squared and integrated over the momentum transfer according to equation 
(26). 

The following cross sections were calculated using both SWF I 
and SWF XI: Li(2s **■ 2p) , Li(2s + 3p) , Na(3s -> 4p) , Mg(3p -> 4s), Ca(4s -► 4p) 
■and K(4s -> 4p). Two computer programs which were written to do the cal- 
culations are given in the Appendix.. The first program calculates the 



total cross section using SWF I as follows: First the a 's given by 

JK 

equation (23) are calculated. The the b *s equation (24) are calcu- 
■i 

lated. Next, the F T s equation (2 1 ) and the A ' s equation (2 r ) are 

nn n nn 

initialized to 1 and 0 respectively. The I's and A^_.'s are ca ^- cu -*- a_ 

ted next. These must be calculated together since they depend upon 

each other.- Then using the calculated N, T s and A , ! s the F ,'s 

' J nj nj 

equation (2) are calculated. Finally the last part of the program, cal- 
culates the scattering amplitude f (K) equation (22) squares it and in- 
tegrates | f (K) ] over the momentum transfer using Simpson’s rule as an 
approximation to the integral. 

The second program calculates the total cross section using 
SWF II. This program calculates the radial integral in addition to 
the integral over the momentum transfer. Thus, the N . T s and the A *s 
are calculated first. Then the radial integration is done and finally 
integral over the momentum transfer is done. Simpson's rule is again 
used for both of these integrations. 

The results of the calculations are given in graphic form on 
the following pages. These results are compared with the calculations 
of Vainshtein, Opykhtin, and Presnyakov (VOP) and with experimental data 
where available. It was found that the difference in using SWF I and 
SWF II for all of the calculated excitations was less than about 5%. 

The data are therefore only plotted for one set of wave functions. 

Agreement with VOP is very good for the Li (2s 2p) transition. 

Both calculations are high x-rith respect to the experimental data given 
in reference ( 9). Agreement with VOP for the Li (2s 3p) transition 



is not as good with VOP obtaining a peak cross section which is more than 

■ t 

3/2 times and at a lower energy than the SWF result, No experimental 
data could be found for this transition. The SWF result for the 
Na(3s ->• 4p) is about an order of magnitude smaller than, the VOP calcula- 
tion up to 5 threshold units and a factor of 3-8 smaller for the higher 
energies. However, the SWF cross section for this transition has an 
oscillation in it. More will be said about this oscillation and attempts 

to understand it (Z. = 1.35 curve) in the conclusions. As can be seen 

4 

from the graph the SWF cross, section is closer to the experimental values 
given in reference (9 ). . No other calculations or experimental data could 
be found for the Mg(3p 4s) and Ca(4s 4p) transitions. Good agreement 
with VOP is obtained for the K(4s 4p) transition. The cross section 
peaks^at roughly the same energy but the VOP peak is about 50% larger 
and does not tail off as fast. 

All of the Modified Born results show the expected behavior 
of' a lower cross-section at the lower energies with the approach to the 
same value as the Born Approximation at the high energies. 

These data are presented in tabular form in Tables I and II on 
pages 27, 28 and 29. 
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TABLE I 

ELECTRON EXCITATION CROSS SECTIONS 

? 

( Units of Tra^ ) 



Li (2s -*■ 2p) 


Li (2s -> 

3p) 


E ' 

(threshold units) 

Born 

(SOT) 

Hod , 
Born 

(7) 

VOP v u 

cc (8) 


Born 

(SWF) 

Mod . 
Born 

VOP ^ 


1.16 

.101(3) 

.961(2) 

.998(2) 

.284(2) 

— 

.551(0) 

.399(0) 

.201(1) 

I 

1.64 

.130(3) 

.125(3) 

.131(3) 

.564(2) 

— 

,788(0) 

.518(0) 

.187(1) 

f 

i 

j 

2.00 

.127(3) 

.122(3) 

.129 (3) a 

— 

.159(2) 

.917(0) 

.667(0) 

•151(l) a 

f 

i 

}. 

2.44 

.119(3) 

.115(3) 

. 122 (3) a 

— 

.230(2) 

.102(1) 

.813(0) 

. 127 (l) a 

\ 

\ 

5 

3.56 

; 

.995(2) 

.966(2) 

.103(3) 

— 

.239(2) ' 

.118(1) 

.101(1) 

.891(0) 


5.00 

.818(2) 

.797(2) 

. 850 (2) a 

— 

.213(2) 

.120(1) 

.107(1) 

.631 (0) a 


6.76 

.674(2) 

.659(2) 

.692 (2) a 

— 

,186(2) 

.115(1) 

.105(1) 

.477 (0) a 


11.24 

.475(2) 

.465(2) 

.501(2) 


.161(2) 

.976(0) 

. 911(0) 

.288(0) 


17.00 

.351(2) 

.345(2) 

.375(2) a 

. 

.156(2) 

.809(0) 

.765(0) 

, 190(0) a 


•: 24.04 

.271(2) 

.266(2) 

. 298(2) a 

* 

.144(2) 

. 

.675(0) 

.643(0) 

. 133(0) a 



a - Graphically interpolated. 




TABLE I 
(cont . ) 


Na (3s -»■ 4p) 

Mg (3s -> 

4p) 

K (4s -v 

4p) 


Ca (4s -> 

4 P ) 

Born 

(SWF) 

Mod. 

Born 

VOP ^ 

Exp^ 

Born 

(SWF) 

Mod . 
Born 

Born 

(SWF) 

Mod. 

Born 

VOP 

Born 

(SWF) 

Mod . 
Born 

.266(0) 

. 211 (-1) 

.320(1) 

— 

. 753(-l) 

. 195 (-1) 

.103(3) 

.417(2) 

.135(3) 

.411(2) 

.294(1) 

.211(0) 

. 187 (-1) 

.360(1) 

.232(0) 

.214(0) 

. 842 (-1) 

.130(3) 

.597(2) 

.187(3) 

.774(2) 

.112(2) 

.194(0) 

. 245 (-1) 

.347(1) 

.234(0) 

.253(0) 

.119(0) 

.125(3) 

.611(2) 

. 186 (3) a 

.866(2) 

.161(2) 

.194(0) 

. 381 (-1) 

.311(1) 

.230(0) 

.273(0) 

.148(0) 

.115(3) 

.596(2) 

. 182 (3) a 

.902(2) 

.205(2) 

.209(0) 

. 240 (-1) 

.222(1) 

.223(0) 

.279(0) 

.178(0) 

.931(2) 

.‘529(2) 

.159(3) 

.873(2) 

.264(2) 

.219(0) 

.103(0) 

.170(1) 

.214(0) 

.261(0) 

.183(0) 

.746(2) 

.452(2) 

•135(3) a 

.786(2) 

.288(2) 

.219(0) 

.120(0) 

.132(1) 

.193(0) 

.236(0) 

.175(0) 

.604(2) 

.383(2) 

. 114 (3) a 

.690(2) 

.289(2) 

.196(0) 

.124(0) 

.839(0) 


.187(0) 

.148(0) 

.414(2) 

.280(2) 

.821(2) 

.526(2) 

.261(2) 

.169(0) 

.120(0) 

.585(0) 

— 

.149(0) 

.122(0) 

.301(2) 

.212(2) 

. 595(2) a 

.409(2) 

.225(2) 

.144(0) 

.109(0) 

.426(0) 

— 

.121(0) 

.102(0) 

.230(2) 

.166(2) 

.457 (2) a 

.326(2) 

.193(2) 


a ~ Graphically interpolated 

b - Additional points calculated, n = 1.02, a ~ 0.148; n - 1.04, o = 0.195; n = 1.03, o = 0.241. 




TABLE II 


ELECTRON EXCITATION CROSS SECTIONS FOR Na (3s -> 4p) USING Z = 1.35. 

2 

(Units of na )-. - 


E (Threshold units) 

1.16 

1.64 

2.00 

2.44 

3.56 

— 

5.00 

6.76 

11.24 • 

17.00 

24.04 



Born (SWF) 

.469(0) 

.477(0) 

.555(0) 

.657(0) 

.818(0) 

.883(0) 

.878(0) 

.781(0) 

.664(0) 

.563(0) 

Mod. Born 

.313(-1) 

. 581(-1) 

.119(0) 

.202(0) 
, - . , 

.367(0) 

-.478(0) 

.532(0) 

.537(0) 

.489(0) 

.433(0) 
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The line strength for spontaneous emission in the dipole 
, 10 

approximation is given by 

n.itL. 2 

S = l |i 1 k | (27) 

i,K 

where i and K number the degenerate sublevels of the upper state n 
and'j lower state m; and where the summation is over all possible 
combinations of the sublevels of the upper state with those of the 
lower state; and where 

J K™ = <n|r|m> (28) 

It will now be shown how S is related to the Einstein coeffi- 
cient for spontaneous emission. This will be done semiclassically in 
the sense that the E and M field is treated classically while the par- 
ticles are treated quantum mechanically. ' A correct quantum electrodyna- 
mic treatment leads to the same result s^. 

Maxwell’s equations of motion for the electromagnetic field 
are in gaussian units 


V * E = 4t rp 

(29) 

-y -y- 

V x H - 

1_ 

c 

SE 4tt 
at. ~ c J 

(30) 

-v 

n v a. ^ ^ n 

V X E + — = 0 

c at 

(31) 

V * H = 

0 


(32) 



31 


From equation (32). we sea H = V x A. If this is substituted 
into equation (31) and the order of the spatial and time derivatives 
is interchanged we obtain 


v x (E +-i ff) = 0. 


Thus 


or 


E + i = - 7 <j) 
c at v 


1 3 A. 

E = -- -- V* 


Substitution of the equations for E and H into equations (29) and (30) 

V(7-A + i|f)-V 2 A + i — = ^ J (33) 

c 3t 


1 3 Tb 't , . 

— — — V • A + V 6 = - 4mp 
c 3t 


(34) 


where V x(V x A) = V (v • A) - has been used. Nov choose the Lorentz 

gauge so that 

c 3t 


Then equations (33) and (34) become 


2-> 

V 2 t-L. LA 

\ A 2 2 

c 3t 


4£ J 

c 


(35) 


n 2 , 1 3 q> / 

V <P j — j = 4 

c 3t 


(36) 
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Now, taking the curl of equation (35) we obtain 


zj 1 8 H 

V H 2 — k 

c at 2 


A t r 

- — V x J 
c 


(37) 


Assume that the three cartesian components of J vary harmonically in 
time with the same angular frequency to but not necessarily with the 
same phase: 

J (?, t) = 2 [j <?)| cos (lot - 8 ) = J x (?) + J* (?) e i “ C 


where 


i0 

J (?) “ l J e 


with similar equations for the y and z components. We are interested 

-y —y 

only in the solutions for E and H that have the same frequency to. Thus 
the x- components are given by 

E^(r , t) = 2|E x (r)| cos (tot - = E^(r) £ 1Wt + E x ( r ) e 1Wt 09 ) 


H (r , t) = 2|H"(r)| cos (wt - a ) * H (?) e 1£i)t + H* (?)e lLJt (40) 


where 


x 


E (?) - 
x 


E (?) e 1 ** 

X 


and H (?) = |H (l e la x 

X 1 X 


Using equation (30) E is given in terms of H in a vacuum by 


E(r) = — V x H (r) 

to 


(33) 


( 41 ) 
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Substituting the expressions for J, H and E into the wave equation for H, 
equation (37), we; obtain, 


(v 2 + K 2 ) H(r) - - “■ 


V x J (r) K = lu/ c (42) 


The retarded solution, of the above equation is given by 

1 f V x J (r 1 ) iK|r - r 1 | d 3' r , 


-4- 

H (r) = 


[-V -> 

r - r 


<43) 


Now 


iKlr - r* 


!»-*- .->• t 

r - r 


e iK(r - r' cos 0) (44) 


where 8 is the angle between r and r* . 

The Poynting vector is 

^ g(?,t) x H(r , t)] 

Using equations (39) and (40) for E and H, we see that the average over 
a period of oscillation has a z component 


(?) = — ^|E x (r)| |H y (r)| jcos(wt - d x > cos (wt - a y ) 


time 

ave 


- E (r) H (r) 1 cos (wt - <f> ) cos (wt 
v ' x y 


- “xil 


time 

ave 


P z (r) 


c 

2ir 


j|E x (?> | |H y (r) 


cos (<j> - a ) - E (r) E (r) cos (a - 4 

Y x y y x x ; 
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Then 

P(r) = Re [e(t) x (45) 

Now we are only interested in those terras in the energy flux 
o 

that fall off as 1/r since only these correspond to radiated energy. 
Choose the z-axis along the vector r" which goes from the center of the 
charge-current distribution to the point at which the field is measured. 
Using equations (41) , (43) , and (44) we obtain 


7T iK lKr / -lKz* d r’ 

H -*■ ■ e I J (r ) e 

x re- J y 


r 3 

„ iK iKr _ -iKz d r T 

H *+■ — e J (r') e 

y rc x 


H 0 (46) 

z 



E' 

y 




E -4-0 
z 


Partial integration has been used to get rid of the derivatives 
of ” ... . iricegran'd of H. . The above equations show that the asymptotic 

fields are transverse to the direction of propagation. They also show 
that only the component of the current perpendicular to the direction of 
propagation contributes to the radiated energy. 
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Substituting these fields into equation (45) we obtain 

r 


p « 
z 


K 


o 2 * 
2irr c 


J e 
x 


iKz 1 ,3 , 


dr' | + | J e 

1 1 y 


-iKz’ ,3 , 


d r* 


Generalizing the above equation, to give the average energy flux in the 
direction of K we obtain 

2 i 


P = 
z 


K 


' 2 
2rr c 


V r,) 


ill* r * 


,3 t 
d r 


where J v is the component of J perpendicular to K. 

l k 

In most cases of practical^ interest the wavelength of the 

radiation is many times greater than the linear dimensions of the cur-* 

-iK* r 1 

rent-charge distribution. This means that R*r' «1 and e can be 

replaced by 1 in the integrand of the equation for P , Then 


P = 


K 


- 2nr 2 c ' -L K 


I f J 1Ir £'> d 3 r' I 


(47) 


From equation (46) with the same approximation it is apparent that the 
polarization of the radiation is determined by the total current vector 


J = 
o 


3(?’) dV 


*7" 

With J q linearly polarized the integral in equation (47) can be replaced 

by ' ^ | sin 6 where 0 is the angle between K and J . Then the total 
J o' o 

power radiated is 

^ 2lT TT 


P = 


K 


27r r c 


5 l*ol 2 



r^ sin^B de d<f> 


o o 


2 2 

|J I 

5 o 1 


(48) 


3c 
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Equation (48) is also valid for non-linear polarizations. 

-*v 

To convert to quantum theory the classical expression for J 
must be replaced by its quantum analog. We want to replace J by a 
current density that is associated with an initial upper state and 
a final lover state u^, since energy is radiated during the transition 
from K to n. It is natural to represent the current density as a pro- 
duct of a charge density and a velocity and to take for the velocity 
the momentum operator divided by the mass: -(ih/m)V. The charge den- 

sity for a stationary state is expected to be the charge of the particle 
times its position probability density i.e., e|b| • However, we are 
concerned here with a transition between states and so replace this by 

- 4 - , 

eu u . Thus we can replace J(r) by 
n K. 


deft * t*\ 
j( r ) + - u (r)Vu (r) 

m n n . 


To get J we must integrate the above expression for J(r) over the coor- 


dinates. 


Thus 


-> _ _ ieffi 

o m 


* -v 3 
u V u„d r 
n K 


or 



* + 
u P 
n 


,3 

U K d r 


or 


J - - <n |P| K> . 

o m. 11 


(49) 
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-V -j- d r 

The matrix element of P is F = m-rr . 

dt 


1 i d “>■ 

So - <n P K> = <n r K> 
m * at 


If an energy representation in the Schrodinger picture is used we have 


<n(t) [r K(t)> - < n(o) r K(o)> e 


iw t 
nK 


E ~ E v 
n K 


, ia) v t 

~r~ <n(t) | r j K(t) > = t: <n(o) |r| K(o)> e n 


= ioo -,<n(o) r K(o)> e 


iw v t 
nh 


= iuj <n(t) |r| K(t)> (51) 


Substituting equation (51) into equation. (50) we obtain 


r <n l^l K> = ±bi <n | r [ K> . 


Substituting equation (52) into equation (49) we obtain 


^ = ieta v <n r K> 

o nK 


Substitution of equation (53) into equation (48) gives the radiated power 
We can interpret this power as the product of the spontaneous rate of 


transition from n to K and the quantum energy dioo v = E v - E given off 

ni\ a. n . i- 

in each transition. 
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Thus 


P « 


4K' 

3c 


„ ! <n I r I K> 

nK 1 1 1 


and the transition probability per unit time for spontaneous emission 
(Einstein probability coefficient) becomes 


, 2 2 
4K e u 


nK 


A - ' 

nk 


<n r [K> 


Using a) ~ Kc we obtain 


A nK " 


3iic 


. 2 3 
4e w 


nK 


3ilc‘ 


<n r K> 


The line strength is the sum over all possible combinations of the 

(possibly) degenerate sublevels with those of the lower state of the 

matrix element in the expression for A Tr . 

nK 

To calculate the line strength of the Na(3p -> 3s) line we 

need the effective atomic numbers Z for sodium. Using Slater’s rules 

n 

we find: Z - 10.7, Z^ = 6.85 and Z^ = 2.20. From these values the A’s 

/ 

and N’s can be found and substituted in Slater’s radial wave functions. 
The same result is obtained for the line strength of the Na(3p -»■ 3s) if 
either SWF I or SWF II are used. 

Thus we obtain (in atomic units) 


for m = 0 


& -y "5 

i(j r if>__ d r - 2.76 e 
n • ‘K z 


for m = 1 


* 3 / — ■ i * A 

\jj r t|» d r = ,799 JfT -e + ie 

n K x y 
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r 


/V A 


for m = 1 


ib r d^r = .799/6* fe + ie 1 . 
v n K L x y 1 


When these results are substituted into the expression for S the answer 
is (in atomic units) 

S - 22.9. 

This compares reasonably well with an experimental value of 19 atomic 
units given in reference 2. 

The value obtained for the Einstein probability coefficient 
is 


^nK = 2.256 x 10^ sec 





CONCLUSIONS 


With thd exception of Na(3s -> 4p) cross section the SWF cross 
sections are roughly with a factor of 3-4 of VOP’s results- In parti- 
cular, there is very good agreement with VOP for the Li (2s -f 3p) cross 
section. One should realize that comparing one theoretical calculation 
with another does not establish the validity of either. 

The oscillation in the Na(3s 4p) cross section is discon- 
certing. Normally such oscillations come from the interference of coupled 
states, but couplings haye not been included in the SWF calculation and 
so it was first thought to be a programming error. A second program was 
written to serve as an independent check on the results of the first pro- 
gram. In this program the radial integration was also done on the com- 
puter. Numerous excitation levels of hydrogen using hydrogen wave func- 
tions were calculated with both programs. Both gave the same results 
and agreed with the results given in VOP. When SWF's were used the pro- 
grams also calculated the A and N parameters. Both programs used the 
same routine for calculating these parameters. So it was concluded that 
the oscillation in the Na(3s 4p) cross section was real. 

By setting all of the A's individually equal to zero it was 
found that only A^.^ significantly affects the cross-section. In -par- 
ticular by setting A^.^ equal to 90% of its calculated value the oscilla- 
tion is removed. However, the cross-sections increase by about 60%. 

Using the orthonormality conditions <R , ,R f > - 5 , it can 

G n£ n Z nn 


be shown that 
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413 


drr 2 R 31 (r) [r 3 exp(-Z^r/4)] 


Thus, the absolute value of can be reduced by using a different 

value for Z^ than that given by Slater's rules. The value obtained 
for by using Slater's rules is = 1.00. By setting Z^ - 1.35 the 
oscillation is removed. However , this increases the cross section to 
3-4 times the original calculated value and the peak shifts toward 
higher energy. This rather larger change in the cross section should 
be expected since theA's, N's and F's all depend on Z^ and these 
appear in the expression for the scattering amplitude f(K). 

Thus it appears that the absolute value of t0 ° large. 

This could be due to too large an overlap between R^ and the. characteris- 
tic part of which is the factor in brackets in the above expression 

for \l3- 



APPENDIX A 


The following programs were used to do most of the calcula- 
tions related to this work. The first program does both the radial . 
integral and the integral over the momentum transfer using Simpson’s 
rule as an approximation to the integrals. The first part of the 
program down to and including statement 10 calculates the N atl ^ 
^n£j' S * ^en ra dial integration is done and finally the integral 
over the momentum transfer is done. 

The second program as given uses the SWF I. Only minor 
modifications are required to use SWF II. In this program only the 
integral over the momentum transfer is done as equation (1) is used 
for f(K). This program will work for any s **- p transition for which 

n,. > n.. If n '< n. enough N.’s* A ,'s and F , *s may not be calculated 
f - i f x r j nj nj 

Only the statements marked * need be changed to calculate the total 
cross-section for other than s p transitions for which > n . * 

t - i 



\ r. 


/ / c 1 •>. “ c: r jo Si ** 1 6 '? ? .< r; - , ^ v - a i j l g i s k ■ , i . p l - 

/ / e: x e c r o r t c- c l g , p a o n = m c m a p ■ 

//FOOT , C .vc. t m jf) * 

0 1 MF N S I ON T N ( 1 0 ) , X SUN < F ) * 4 X ( 3 ) , A ( 3 ) , F >: { S') , X A N 5 ( 3 ) , 5 

D I Y F NS I ON P F X ( 5 ) , R I F A ( 2 1 C ) » R ! < 2 0 3 ) » f: ’ F ( 2 0 3 ) 

0 ] NON S J ON C ( 4 * 4 ) , A A ( A * 4 • A ) * R ( 4 , 4 * 20 3 ) * ? { 4 ) . 

0 f; A 0 3 0 , NT , NF , { 1 ( I) ? 1 = 1 , K F ) 

SO FORMAT (?!?» 10F3.2) 

•OF AJ( 5,50) ( A ( 1 ) , r~ 1 * 5 ) , ( T N f I ) > I =1 *3.0 

5 0 • FOP M A T { A ( F S „ 2 , 2 X ) / 1 0 ( F 6 . 2 , 2 X ) ) 

? A' FOR U A T ( F ]. A . B , 4 X ,11* A- X * 1 1 ) 

A i, form AT ( F 3 6 * B , A x , H 3 , 4X % I 3 ) 

7 7 ~ 1 1 « ~ r 
T i: =e>7U'0 
A A A - 1 * ' 7 7 ? 

CO -6 «>-•- A. A A * ? / T T 
XN - 3 0 0 * 

NX " X K + 1 c 
F?N = 200. 

N P “ \< N •«- 1 * 

00 3. I = 3 * N F 

00 1. J = 3,N C 

no i K - i * n f ; • 

A A t T * J ,Jn ) -* 0 » 0 00 i 

1 CONTINUE: 4 1 

DO 3.o JO , nr 

no 3 o l •• 1 j N r 1 

IF (Ji-l-L) 3 0*10,30 ' 

IS A S 0 0 ~ J o Ei 

1 F C J ■■■• 2 ) 4,2,2 


) 


£- 


2 

A 

4 

5 

6 
7 

p 

a 

10 

ADO 


K = J- 1 

DO 3 1=3. 5 K 

A.SOR = ASOR + A A (J,L, ! > **2 
5" 2'!* J + 1 
G~ J 

C( J •, L } = 1 * / ( SORT < GAMMA ( S ) * ( G/ ( 2 . *2 ( J ) 1 ) **{ 2*J + 1 ) -A: 
WRITE (6*24) C(J>1> » J>L 
DO <? N ~ 2 ? OF 
H- C N-J ) <=> , o , *3 . 

ASUM=C.O 
IF (J-2) S t 6»5 
l. L" J~ 1 
DO I-1,LL 

ASUM--ASUM+AA C N * L , T ) *AA 
SA“M*i-J + l 
ft — J 
Q = N 

AA ( M * L , J ) =C ( J , L ) * ( A$UM“ GAMMA ( SA ) * ( U / ( l { J ) /H+ 7 ( N ) /( 
AA ( /( t 2 . 3 ) " . 96*66 (4*2*3) 

V'R I T £ (6,34) AA ( N * L , J > , N . L, J 

CONTINUE 

CONTI MU ! : 

.PRINT 4 00 
FORMA! (1H1) 

A 1. 1 UP - 4 

0 0 L P = ( DL r-’R-ALl MR ) / RN 
V ~ A L T'0 P - D R L 9 

M") 60 I = 1 j 2 ^ j ' 

DO 6 0 IF1FF 


■ OR ) ) 


) ) * ( N J 1 ) ) 


.c- 

4 > 





no 

60 

L 

= } 

F! :: 












60 


R ( N 

9 0 9 

l 

) - 

r. 

' r « s' 

JO 













R( 1 

) 1 ) 

l 

) - 

C( 1 

,1 I 













Ir 

< NF 

- 

2 ) 

6 3 

s 62 


6 2 









6 2 


R ( 2 

•; i , 


) « 

C ( 2 

, 1 . ) 

>; 

AA( 2 9 1 

9 

1 >* 

C ( X 

9 i ) 







IF’ 

{ N " 

- 

7 ] 

6 3 

9 6 7 

• 

6 3 
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R ( 3 

) i i 

1 

) 

C ( 3 

, 1 ) 

* 

{ A A ( 39 

1 

92} 

*C < 

2 , 1 ! 

I ' ,r HA ( 2 9 1 9 1 ) 

*C ( 1 

9 1 ) + A l\ 

( 3 9 1 9 1 } *C ( 1 9 1 


1 . 

1 ) 


















IF 

( NF 

- 

4 ) 

6 6 

. 6 

9 

64' 
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R ( 4 

> 1 , 

1 

i. 

) “ 

C { 4 

? 1 I 


( A A ( 4 , 

] 

9 3 ) 

*C( 

39 1 ! 

1 * A A ( 3 • 1 9 2 } 

•' C ( 2 


( 2 9 1 9 1 ) * 0 I J 9 1 


1 

) + A 

A ( 4 

i 

1 , 

2 ) * 

C( 2 

♦/ 

3 J *AA ( 

2 

9 1 5 

1 ) * 

C ( 1 - 

! 3. ) + A A ( 4 9 1 9 

1 ) *c 

( 1 9 1 ) ) 

. 


6 f> Y •- Y -i- 0 '£ L F > 


DO 7 0 J " 2 5 M 7 

v' -• Y -\- p F l R 

DO e-0 L “ 1,2 
no B 0 N - 1 , R F 

f.- 1 -j .. f,t _ i 



IF ( J- 3 6 j 8 2 , 4 0 

,40 

4 0 

I F ( N-l 1 80 ,30 

ft 'i 

6 3 

IF (J~-; r) 5 R 2 j 4 

1 9 4 1 

41 

] F ( M - 2 ) BO t 8 ^ , 

O 9 

o 2 

DO 61 K = L * "11 


81 

H ( N 9 L 9 J ) - • R ( N 9 L 9 

J ) + C ( N 9 L ) * A A ( K , L , K) * R ( K 9 L 5 J ) 


IF (N-l) 

V '.7 : 

9 4 

R ( M , L 9 J } ” C ( N L } 

j- F X R { ~Z ( N J >-• Y / R ) + R ( A 9 L 9 J ) 


GO TO 60 


9 3 

R (N. 1 . 9 J ) = C ( N 9 L ) 

* ( Y * * { H ~ 1 ) -x- R ;< p ( - ;/ { N ) y / N ) ) -{• R. 

8 A 

COOT iiMJJ 



R I ( J } " R ( 3 , 1 , J ) 


7 0 

Rr( J) 4 ,2, J) 
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